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ABSTRACT 


A theoretical and experimental investigation of the stability of 
developing laminar flow in a circular pipe subjected to small, rota- 
tionally symmetric (torsional) disturbances is presented. Analytically, 
developing mean velocity profiles were represented in polynomial series 
form and subsequent formulations yielded determinations on stability 
which reduce to the classical conclusions for the case of fully developed 
flow. 

Experimentally, torsional disturbances were imposed upon a developing 
laminar pipe flow by sinusoidally rotating a 1.5 inch (1 diameter) length 
of pipe in the entrance region of a carefully constructed wind tunnel 
facility. For flows with maximum Reynolds numbers of 8500, maximum 
amplitude and frequency of oscillation were 4% inch and 30 Hz., respec- 
tively. Hot wire anemometer measurements verified the stability of the 
developing flow field. Moreover, measured wave speeds and decay factors 


compare favorably to caiculated values. 
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NOMENCLATURE 
Constants, coefficients of the power series 
solution of the disturbance amplitude function 


Constants, coefficients of the polynomial 
approximation of the velocity profile 


Dimensionless wave speed 
Amplification factor 

Inside diameter of the pipe 
Pressure 

Radius of the pipe 

Reynolds number (Rows) 

Mean velocity in the pipe 
Maximum velocity in the pipe 


Dimensionless velocity profile in the entry 
region, referred to ve 


ax 

Axial component of velocity 
Tangential component of velocity 
Axial coordinate 

Disturbance wave number 


Kinematic viscosity 


Stream function 
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I. INTRODUCTION 


The purpose of the work reported here was the study of the evolution 
of small, rotationally symmetric disturbances introduced into the laminar 
transition region (the hydrodynamic entrance length) of Poiseuille pipe 
flow. 

Since the pioneering work in the field of hydrodynamic stability 
conducted by 0. Reynolds in the 1880's, the stability of fully developed 
pipe flow has been the subject of many theoretical and experimental 
studies. All of the analyses up to the present concerning the stability 
of fully developed pipe flow show that such a flow is stable to small 
disturbances. 

A review of the theoretical and experimental determinations con- 
cerning the transition phenomenon leads one to the following considera- 
tions: (1) transition impends because the disturbances introduced 
contain some initially large disturbances so that small disturbances 
are not significant regardless of when or where they are introduced; 
or, (2) initially small disturbances, either axially-symmetric, rota- 
tionally symmetric or three-dimensional, grow in size along a favorable 
length of conduit, namely the entrance length, via convective inter- 
action while the boundary layer is rapidly developing. Consequently, 
these small disturbances, which would not have grown in a fully developed 
region, now acquire a magnitude which eventually renders the fully de- 
veloped flow unstable. 

A theoretical analysis in support of the above considerations has 
been presented by Tatsumi [1]. His conclusions demonstrated that small 


disturbances magnified in the entrance region of a pipe flow field above 
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a minimum Reynolds number of 19,400, do cause transition after a passage 
of 17 tube radii. The significance of his study is not the determination 
of a value for the critical Reynolds number, but that there is a pre- 
diction of instability. 

Presently, the case of developing axially-symmetric flow is treated. 
Results are shown to be reducible to the determinations of Pekeris [5] 
for fully developed pipe flow. Additionally, experimental verification 


of calculated wave speeds and amplification factors is presented. 


iz 


II. THEORETICAL ANALYSIS 


Ae CYLINDRICAL ORR-SOMMERFELD EQUATIONS 

To describe a general disturbance in mathematical form requires the 
formulation of a three-dimensional (r, @, z) vector function of time, t. 
The generally accepted treatment for small disturbances is to view them 
as a Fourier series of a large number of elementary disturbances, each 
of which can be represented in time and space by a finite number of para- 
meters. Indeed, this is the foundation of the small disturbance theory - 
that the disturbance is a linear combination of spatially periodic 
components which grow or decay with time. Such a component is character- 


ized by the stream function, 


Y= F(r) Exp a + (Gt —Cx2| (1) 


By further contending that the imposed disturbance be rotationally sym- 
metric, then the @ dependence is eliminated, and the representation 


reduces to, 


Yo = Ftv) Exp |Cx(ct-2) (2) 


Therein ©C is a real number - the wave number of the disturbance in the 
fluid. C is a complex number, c* ae iC. . The sign of the imaginary part 
OFC, Ci» determines whether amplification (Cc. < 0) or decay KC 0) 
occurs, whereas C. is the wave speed of the disturbance. Experimentally, 
such a disturbance is difficult to generate, since, as formulated, it 
should be imposed upon the entire length of the conduit over the entire 
cross-section of the conduit simultaneously, and very periodically in z 


as well. 
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Analytically, the disturbance described in (2) is thought of as 
being imposed upon a steady mean flow characterized by a stream func- 


tion, 
4. & Yl Ha) (3) 


The equations of motion for incompressible flow symmetric about the 


axis of a pipe are: 


Br ve uae = LRP Pe] ww 


Sey Ber UM m5 +7FU ia, 
SAm T AX - 6 (4c) 
a am ot 7 

ao = >, 
v= = oo 2. ae +=, =} = (4d) 


5 2% 


It is convenient to introduce the non-dimensional = 
(a coh siy= £U 
ease » t= ie? = ie 

and define R, the Reynolds number as R = us, With these substituted, 

the system of equations (4) becomes: 


DA el OM RU on oe +4 i 
Sut vw” Ses vw - ei 


/ Md 2 
v5 - +4 View) 
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| a 
a te Sivan aun (Sc) 


Defining U = 13} and V = -t2¥ » it is noted that with the dependence 
on @ dropped, that there are two possible types of symmetrical distur- 
bances which may be generated and which are not mathematically coupled. 
First, there exist the torsional disturbances, which are concentric 

with respect to the pipe axis; secondly, there exist the meridional 
disturbances, which occur in planes passing through, and parallel to, 

the axis of the pipe. Treatment of each case is necessary due to the 


aforementioned difficulty in experimentally imposing the mathematically 


exact disturbance upon a flow regime. It is assumed that, 


V(r 24) = f(r) Exp|éx (Ct -2)| (6a) 


N(wzBet) = GVEXp[éa (ce-2) (6b) 


These relations, together with the definitions of velocity components 
in terms of the stream functions, lead to the following uncoupled sys- 


tem of exact equations of motion: 

s pe 82 oY & Hy 7 

SOY +a so+ they 33? ¥ ay ay] 4 (7a) 
SOY = 4 DIY 

2 i ed ow 3 = 

BL 4 Lae a4 - ai £ per (7b) 


where a AC fs >. 23 
a eS 


As Y is composed of the main flow with the disturbance superimposed 


upon it, it can be described as, 


‘oe Yo + Y’ (8) 


/ 
where Y, refers to the main flow and Y is as defined in (6a). 
Substituting (8) into the system (7), dropping squares and products 
of yw! and {2 , as each is assumed small, leads to the following 


uncoupled system of linearized equations: 


Sey + Ef Bee SO yy BY ROH - Se Pry aE 3, Dillon 


Pe Oye Deep = ee 


R. 
> aY atk _ o Ot 4 pr 
a5 + c yr 82 ce oY R DIL, (9) 


V\ 


To proceed further requires that the exact nature of A be known. 

In the fully developed case, Y is a simple function of ff alone, 
and is easily treated. In the case of developing flows, such as the 
one under investigation, there is no precise, known functional descrip- 


¥, 


ELOnwOL. 


l\ 


Yo (132,4) | 


Ds DEVELOPING FLOW STABILITY DETERMINATIONS 

Schiller [2], and Langhaar [3] present workable models of developing 
velocity profiles in parametric form. However, after long and fruitless 
attempts, the author has concluded that there is no simplified mathematical 
solution suitable for representation of the developing mean flow at pre- 
sent. Accordingly, following Tatsumi [1], purely parallel flows which 
have the same velocity profile at a given axial distance as the developing 
flow shall be treated. For such velocity profiles, it was found that 


experimental data could be approximated to an excellent degree of 
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meee ee 


accuracy by a third degree polynomial obtained from a least-squares 
curve fitting procedure. To make the solution more general, it is 


assumed that 


U, (r) be thar thr +b, rr? 0 


from which by 
= 3 , > 
Yo lr) = boyy barry be 4 pba br > ay 
2. > 4 Ss 
is obtained. By obtaining the necessary terms from (11) and (6) and 
substituting them into (9), the following equations are obtained: 
a = 
D*g = CXR[( Cob, - bir-b,P—bsr*p H3h- ai (12a) 
- 
where as | f! fe (12b) 
J = f —- x f 
and primes denote differentiation with respect to az Treating (12a) 
first, it is clear that for large R, the terms on the right side of the 
equation will predominate, since even when & is assumed small1,XR is 
> 
still large. Yet, if the entire D4 term is ignored, it is impossible 
to satisfy all of the’ boundary conditions of the problem. Order of 


magnitude considerations allow for the retention of the g term, and 


thus satisfaction of all boundary conditions. Returning to the defini- 


tion of g, 
g's fT Leh s (Rwy af 
Inserting (13) into (12b) then yields, 
gi 1 §"+(%- ar. R %)F"+(- 4S, +Ba\e (14) 
tfORG-cxRF) £ =O 


where 


Og = Ga ( C-bo~b, FH bzb=b, -) 


i (13) 


and 


Sg = La (Sher - bi) 


r 
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Before simplification, equation (14) is of the general form, 
}Vv it] 1 
Yor Pr) go + Qir) gl + RW) Y + Tor) 4 =O 
a form which possesses a regular singular at the origin, r = 0. This 
fact guarantees the existence of at least one series solution of the 
cY\ = K 
form 15 
fir)= 7 Boer = 


Krod b> = BR 
which will converge for0Cr< ro where ty is the distance from the 


origin to the next nearest singular point. As the disturbances intro- 
duced are confined to the vicinity of the pipe wall, r “1, and again 


following Schlichting [6], the terms 


yl }/ 2 il / 
6 f 5 ty ae 5 ~ wf and +, 7 


are neglected as they are all of the order 1, whereas the oR term is 


assumed quite large. The resultant equation is 


pie Reags R 5 + (X*RE-KE\E 008 


Assuming a solution of the form (15), differentiating and inserting 


terms, 


= = 
2 brit Kk om ter! im ee-2 meie3) Ag OR? 
t i va 
— lak ; os lt KK mtie-] a i ia (C=hy- br-by=byr) 
~= Go dine 


pcoek [ Seren) ae peel C-be~b, rb her”) 
Kzo 


+0%«R f = le (<0}(Cbe-h r-bir—byr?) 
bed 


~ (Bb3r-- 2) =O 
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Setting the coefficient of the highest power of r equal to zero results 


in the indicial equation, 


Mm(m-41 \n~-2\m-3 )% = 0 


with indicies n = 0,1,2,3 to guarantee that an 0. 

Before proceeding further, specific attention is due to the boundary 
conditions of the problem. The no slip condition at the wall guarantees 
that f'(1) = 0, tien, + 3450. Since the flow is axially symmetric, 
v' must be symmetric, hence zero on the axis, or £(0) = 0. Finally, 
there can be no sources at the wall, and therefore Wi o1) = 0, or f(1) 
= 0. Here, it should be noted that most previous works also include 
a condition that u'(0) = 0, or that £'(0) = 0. However, as noted by 
Leite [4], there is no guarantee that only those disturbances generated 
by the experimenter are being accepted by the flow, i.e., there is no a 
priori reason to expect u'(0) = 0. It is noteworthy that it is not neces- 
sary to impose this artificial condition in this analysis. 

Taking the largest exponent of r first, it is found that the second 


highest ordered coefficient, 


(2-1 Ymjrrrd \rrre )% =O a 


can only be identically satisfied of a, = O. Continuing, for each co- 


efficient in turn, it is found that the recurrance relation yields: 


~ Cx A (C-hy)(m-2) 
hn (37 tA\(774.2) (™~-2) i 


_ _ ¢xRb - 4)* 
2+ m+ 3\r+rxYne/ A sad 


rs 


= CXR Go [etl ba) ania) _ ralen- iN, 4 o¢'t-by) 
1 Creede Yor!) a7) (7 1) “> 


After substituting n = 3, then; 


GQ, = G2 C&XR(C-bo) 5 Aya —GecaRs, 
i ' r 
ay = Se (rake) _~ L&R (BhrrarC be} 


With these coefficients, the final form of the solution becomes, 


fir) = Ox y> [1+ fs ( xR (C-bo) — eC xR hy (22) 


a ys C+, 
e ee bo) —CxK (3b, + xt >} 


The existence of additional linearly independent solutions is deter- 
mined by returning to the indicial equation. The value n = 0 will yield 
such a solution, but as the boundary condition £(0) = 0 would not be 
satisfied without a QO, then this solution may be discarded. The same 
consideration applies for the value n = 1, i.e., the £'(1) = 0 boundary 
condition would not be satisfied without ae do O. In each of these cases, 
if a, were to be equal to zero, then it is evident that the series could 


be represented as 


R 


n+K+4 
eta 1 
0 
and n = 1 would not then satisfy the resultant 


iY | 


© 


Therein, the values n = 
indicial equation. 
Finally, the case of n = 2, the last of the four index values is 
considered. Now, n = 2 identically satisfies the recurrance relation 
involving ay» i.e., the second highest ordered term. Thus, it is 


possible that this value would yield a linearly independent solution. 
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However, as this exponent differs from the highest by a positive integer 
(vis. 1), it is always possible to generate another linearly independent 


solution of the Bea, 


ot ao 

™ + | mt i 
> KY “+ Mmir) > & 
kK=0 2 


This solution clearly fails the £(0) = O boundary condition, hence it is 
also discarded. Therefore, of the four linearly independent solutions 
possibly generated by the exponents from (18), only one remains that 
satisfies all the conditions of the problem, viz., the solution represented 
in relation (22). 

Applying first the wall boundary condition to this solution yields, 
Cit — eR Ci p SER C: — aoe (6G py c* | 

+ GR canbe ~ db, — cahbr ~ ak (Ce-bo) (29) 

Be a Craig) ce. 0) 


 hEYO 
Separating this equation into real and imaginary parts gives 





Ce~be ~ 2 co aX(p-bo) 3 oie Cr Ce-by } —~O (24a) 
0 


FYO 4540 


ox ! SDP, ee 4 — 
t= am : eat : F590 | Cebt- 08 ] me 


Solving (24a) for C, results in 





SS .-a  aa \aeee 
Ce SR ails) ik ess) SRR 


Sylidem ean F. B., Operational Calculus for Applications, p. 137, 
Prentice-Hall, 1962. 
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Following the same procedure as above in applying the boundary condition 
£'(1) = 0 to the solution in relation (22) yields two additional equa- 


tions 


Ce-bo — fv — bon _ %* (Combe) _ ARC (Ce-be) am 5 (26a) 
5 s 4 /20 /60 


Ree 4 PRE 
= ARS 4 xl (c 2b -G | [=e (26b) 


/2X0 


Solving (26a) for C. gives 


Cran 3O ~ 32h _ 4b: x" a 
KK = ARCe-by =~ KRCg-by 2%R 


It should be noted that C. and - have each been normalized with respect 
to hes and that © has been normalized with respect to d. It is clear 
that the normalized = value, the wave speed, takes on values 0 € CaS 1 
with the maximum value occurring in the region of maximum flow velocity. 
Inspection of (25) and (27) reveals that, within the developing region, 
C, must always be positive. Hence the developing region is stable. 
Further, in the limiting case, viz., in the fully developed flow regime, 
the coefficients of the polynomial representation of the velocity pro- 


file become; 


Analagous considerations demonstrate then that in the limiting case of 


fully developed flow, the solution represented by (22) also nellds. 


i SGe Also reference [5]. 
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Data obtained from the curve fitting procedure is presented later 
to verify the applicability of this solution to the developing flow 
regime. The solution of (24b) would follow from Pekeris [5] under the 
assumption of purely parallel flows. Further pursuit of the exact solu- 
tion of (24b) however, depends on the exact determination of v as a 


function of r and Zz. 
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IfI. EXPERIMENTAL ANALYSIS 


cs TUNNEL FACILITIES AND EXPERIMENTAL PROCEDURE 

The experiments were conducted in a horizontal Lucite pipe 1.5 
inches I.D. and 45 feet (300 diameters) long. A centrigugal fan was 
used to pull air from the stagnation conditions in the plenum chamber 


through the system. The inlet chamber consisted of a rectangular box 





3 feet square in cross-section and 4 feet long which was fitted internal- 
ly with 6-24 mesh wire screens, followed by a 2 foot long settling cham- 
ber. Sizing of the chamber was fixed to limit the maximum interior 
velocity to 1.5 fps. The well rounded inlet horn was a quarter-section 
of a 9 inch by 18 inch ellipse molded of fiber glass and polielee smooth. 
All joints were hand fitted to insure smoothness and uniformity of mating 
interior surfaces. 

The disturbance generator, Figure 1, was located 12 inches (8 dia- 
meters) from the pipe entrance, and was followed by measuring ports at 
16, 20, 32, 40, 48, 56, and 272 diameters downstream of the entrance. 
Uniformity of the internal diameter and concentricity of the pipe were 
found to be within .005 inch for all sections of the pipe. The indivi- 
dual pipe sections were joined by flanges, each pair hand fitted and 
joined to insure uniformity. The pipe sections joined to the distur- 
bance generator were fitted with rubber O-rings to assist in damping any 
vibrations which might carry over from the driver to the pipe. The pipe 
was fastened to a horizontal metal track which was bolted to the con- 
crete deck. The vibrator was securely fastened to the deck after being 
isolated with one inch of corrugated rubber vibration mounting. The 


entire apparatus was longitudinally aligned with a high power transit 
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and leveled. Maximum deviation from true horizontal and longitudinal 
alignment is felt to be no more than .25 inch over the entire length. 

Measurements were taken with a Thermo-Systems constant temperature 
hot wire anemometer utilizing a boundary-layer type probe with a .095 inch 
long, .OO01L5 inch diameter tungsten element. The probe could be located at 
any of the above mentioned stations, and accurately positioned to within 
.050 inch of the pipe wall. 

Prior to the commencement of data-taking, a careful survey was taken 
to determine the physical limitations of the apparatus. It was found 
that the maximum attainable Reynolds number was 9650; however, at this 
setting, the sensitivity of the flow regime was such that any disturbance 
within the laboratory was sufficient to cause transition of the flow. 
The maximum usable Reynolds number was found to be 8500. The survey fur- 
ther revealed that above a frequency of oscillation of 30 Hz., too much 
vibration was being carried over to the conduit. Hence data runs were 
limited to oscillations at 30 Hz. and less. 

The following procedure was used in taking data during the experi- 
mental phase: 

l. The mean velocity corresponding to the desired Reynolds 
number was measured and monitored at the station located in the fully 
developed flow region (272 diameters from the entrance) by regulating 
the amount of air pulled by the fan. 

De The probe was positioned at the station nearest the genera- 
tor, and rotated until the position fixing the maximum reading was found. 

on A velocity profile was first obtained; data taking runs 
then commenced with oscillation frequencies of 30 Hz., with 4% and %& inch 
amplitudes, oscillations then at 20 Hz., with % and % inch amplitudes, 


and finally at 10 Hz at % inch amplitude. 


a 


4. The probe was then moved to the next location down stream, 
and steps 2 and 3 repeated. 

a Upon completion of recordings at all stations in the entry 
length region at the desired Reynolds number, step 1 followed by 2, 3 and 
4 were again repeated for another Reynolds number. 

The following procedures were used in recording data: 

ie The output of the hot wire was displayed on an RMS meter 
and on a dual beam oscilloscope, as was the signal input to the vibrator. 

Dae As the probe was moved from its closest proximity to the 
wall towards the center of the tube in .025 inch increments, the hot 
wire bridge voltage and RMS voltage were recorded and the oscillogram 
either described or photographed. This step was repeated at each loca- 
tion for each data run conducted in (3) above. 

The data obtained in the manner just described was normalized 
through the use of appropriate values. Graphical presentations of this 


data are discussed in the following sections of this report. 


B. TUNNEL PERFORMANCE AND EXPERIMENTAL RESULTS 
Ls Developing Velocity Profiles 

Figures 2 - 5 show the velocity profiles obtained both at high 
and low Reynolds numbers at two points in the developing flow. They are 
compared with the analytical approximations of Langhaar [2], with a cor- 
rection factor of .8 applied to the development length. The correction 
factor was introduced because, similarly to Leite [3], it was found that 
the peak velocity obtained at a given axial location was always higher 
than that predicted by theory. It is additionally noted that the boundary 
layer itself does not develop as rapidly as predicted, although this 


factor is strongly dependent upon Reynolds number and thus may be a 
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characteristic of the apparatus. Figures 22 - 25 show the polynomial 
fitted by a least-squares technique, the results of which were used in 
the determination of the damping factor. The second curve in each of 
these figures is that of a fully developed profile offered for compari- 
son. Table 1 contains a listing of the coefficients of the polynomials 
used. 

The noise level in the developing flow is very much greater, by 
approximately an order of magnitude, than that previously reported in 
fully developed flow experiments. It is felt that due to the type of 
probe utilized and because of it's orientation, some of the noise signal 
is indeed a measure of the radial component of velocity, v. In any case, 
the average turbulence level for all stations was .8 - 1.0%. As might 
be expected, the higher values occurred at close proximity to the wall, 
and the minimum value was found in the core region of the profile. 

The very nature of the hydrodynamic entry region, where most of the 
frictional losses are known to occur in pipe flow, suggest that one 
should expect a rather noisy signal. 

Die The Disturbance 

At the highest frequency and highest amplitude, the turbulence 
level jumped to a peak of 12% near the wall, with the peak at the center 
line of 8% at the farthest station located in the entry length. At the 
highest Reynolds number tested, the level at the station in the fully 
developed region reached a peak value of 1.0% compared with .8% without 
any disturbance. It was noted that at all Reynolds numbers and at fre- 
quencies and amplitudes accepted by the flow, that the radial propaga- 
tion of the disturbance was strongly dependent upon the velocity profile, 
i.e., upon the boundary layer thickness. Penetration into the core region 


was strongly amplitude dependent. 
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a. Reynolds Number Dependence 

One would expect that, for a fixed frequency and amplitude, 
the level of turbulence at any downstream position along the pipe would 
increase with increasing Reynolds number. Figures 6 - 9 show plots of 
Turbulence Level vs. r/a at near and far stations at high and low 
Reynolds numbers and demonstrate this relation quite well. Figure 10 
shows the confining effect of the thinner boundary layer found at higher 
Reynolds numbers, i.e., the higher level disturbances are not allowed 
to diffuse towards the center into the core region. This effect de- 
creases dramatically with decreasing Reynolds number. Not only is diffus- 
ion allowed, but it occurs sooner, as the lower mean velocity would 
indicate. Lastly it is noted that there were combinations of frequencies 
and amplitudes which produced extremely small fluctuations in the mean 
flow, just above the background noise level. Such phenomenon were 
observed at 10 Hz. at an amplitude of % inch and at 20 Hz. at % inch at 
all Reynolds numbers. It appears that these disturbances were simply 
damped within the boundary layer without producing the characteristic 
turbulent spike near the wall. 

Des Frequency Dependence 

Figures 11 - 14 show the very pronounced dependence of 
turbulence level upon the frequency of the disturbance generator. It is 
noted that as frequency is lowered, the speed at which the wave travels 
down the pipe should decrease. Hence the station at which a particular 
level of turbulence occurs for a given amplitude of disturbance at a 
given Reynolds number should be located further up stream as disturbance 
frequency is lowered. This is clearly shown in Figure 15. The idea of 


wave speed varying directly with the excitation frequency would seem to 
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suggest that, since the low frequency disturbance spends more time within 
the boundary layer, it interacts for a longer period with the strong 
viscous forces acting there. Hence it is damped more quickly and com- 
pletely. Again the observation that the velocity profile itself strongly 
effects the characteristics of the disturbance in the fluid is apparent. 

com Amplitude Dependence 

Figures 16 - 19 show the increase in turbulence with 
amplitude at low and high Reynolds numbers. Figure 20 shows the strong 
dependence of diffusion upon amplitude, i.e., at higher amplitudes, for 
a given Reynolds number and frequency of disturbance, the higher the 
amplitude of disturbance, the greater the penetration of the disturbance 
into the core. Lastly, the fact that damping occurs at the frequencies 
and amplitudes tested suggests that there is an amplitude below which all 
traces of the disturbance would damp over a distance cf a few diameters. 
Such a determination was not made here, but its existance is foreshadowed. 
S', The Damping Factor 

Following Leite [4], it is seen that from the definition of the 

u component of velocity and from (6a), 
Uy, me YP (= (x CY oi (28) 
U,, fy 

If a suitable starting point, zo is chosen, this relation becomes, 


uy’ - a C de (29) 
ip EXP L Svs ae] 
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' t(2) 
2.3 LOG, 7 = ~{ ov C dt (30) 
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Thus, 


«; =— 2 Se a 


U! 
OC. dz LOG 10 a (31) 
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If values of LOG, . 7 are plotted vs. the downstream displacement 
oO 


(z - z /d, the slopes of these curves will be, 


4 | 106. {7 
de | HS 
Then, if either a value of wave speed, Co or a value of wave number,  , 


can be determined, then C. can easily be computed. As a check, equations 


25 and 26 can then be solved to yield comparative values. Figure 26 


u! 


log ve (z - z did for data taken at Reynolds 


number of 7600 at the maximum frequency and two amplitudes. As one 


shows a plot of LOG 


would expect, the slopes increase as the center of the pipe is approached. 
Further, at an r/a of .9, the larger disturbance is damped at a faster 
rate. 

Equations 24 were solved utilizing several experimentally obtained 
velocity profiles. Results are presented in Figure 27, which portrays 
the characteristic "hairpin" neutral stability curve. Plots are also 


presented for C, values of .05 and .1 which verify the stability of the 


i 


flows tested. 
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IV. CONCLUSIONS 


(1) Flow in the laminar transition region of a pipe of circular 
cross-section is stable to small, rotationally symmetric (torsional) 
disturbances. 

(2) There exists a strong dependence of amplification factor and 
wave propagation speed upon the shape of the developing velocity profile. 

(3) There is a dependence of both wave speed, Cs and wave number, 


OQ , upon radial position, though evaluation of this relationship will 


J 
require an accurate description of v(r). 

(4) Within the range of frequencies and amplitudes tested, there 
existed a peak turbulence level in excess of 10% very close to the wall 
of the pipe. However, disturbances of this magnitude never reached the 
core of the flow and transition of complete turbulence never occurred. 

(5) Low frequency disturbances are completely damped within the 
developing boundary layer. 

(6) The experimental results suggest the existence of a "critical 
layer'' wherein the disturbance is sufficiently strengthened to cause 
transition to impend. Without this sufficient amplification, the dis- 
turbance diffuses into the central flow region, while continually being 
damped. 

(7) Small disturbance theory per se plays only a small role in the 
overall description of the transition of pipe flows from the laminar to 
the turbulent state. Better understanding of this transition phenomenon 
can only be achieved through an analysis which would describe the spatial 


growth of a finite disturbance rather than the temporal growth or decay 


of an infinitesimal disturbance. However, even the latter, if amplified 


S| 


by the flow regime, quickly reaches an amplitude beyond which lineariza- 
tion of the equations of motion is no longer possible. Hence, there is 
at present no mathematical method suitable for a truely comprehensive 
study of transition. It appears then that further experimental studies 


are the methods whereby more information can be gathered on this subject. 
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V. COMMENTS AND RECOMMENDATIONS 


The experimental apparatus proved functional, easily maintained and 
operated, and very reliable. It is recommended that further pursuit of 
this experiment be preceeded by the following alterations: 

(1) Construction of a new entrance section in accordance with 
the design of Cohen and Ritchie [8] to facilitate flow developing at very 
. high Reynolds numbers. 

(2) Extension of the pipe to allow fully developed flow to be 
achieved at Reynolds number of order 15,000. 

(3) Integration of the air supply into a totally closed system 
to remove the possibility of interference from air currents and tempera- 
ture variations within the laboratory. 

(4) Modification of the disturbance generator itself for use 
of frequencies above 30 Hz. and amplitudes of %'' from the present mechani- 
cal design to an electromagnetic design to eliminate vibration problems. 

After improvements are made to the apparatus, prior to the resump- 
tion of experiments in stability, it is strongly recommended that an 
experimental study of developing pipe flow itself be undertaken to deter- 
mine: 

(1) The developing velocity components u and v as functions of 
the pertinent geometric parameters r and Zz. 

(2) An analytical or numerical solution of the pertinent Orr- 
Sommerfeld equation. 

It is the author's belief that the study of the stability problem 
itself utilizing the small disturbance theory will ultimately lead to a 


better understanding of the problem of beginning of transition from 
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laminar to turbulent flow. Such a study, however, must utilize analyti- 
cal approaches more compatible with experimental procedure than those 


presented in the past. 
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TABLE 1 


Coefficients of the polynomial approximations to the velocity profiles at 
z/d = 16 and z/d = 56 


Reynolds 
z/d number b b b b 


if OZ, = Oe 2 flo -2.63370 


. 733547 -.657861 | 2.64586 ~2.7/294 


823472 -.378369 |1.89644 | -2.34190 


-87101 -~.601644 {2.70048 -3.02695 
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Figure 2: Developing Velocity Profile; z/d 16, R 5700 
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Figure 3: Developing Velocity Profile; z/d 56, R 5700 
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Figure 4: Developing Velocity Profile; z/d 16, R 8500 
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Figure 5: Developing Velocity Profile; z/d 56, R 8500 
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Figure 6: Radial Distribution of Disturbance Amplitude: 
z/d 16, R 5700 
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Figure 7: Radial Distribution of Disturbance Amplitude; 
z/a 48, R 5700 
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Radial Distribution of Disturbance Amplitude; 
z/a 16, R 8500 
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Figure 9: Radial Distribution of Disturbance Amplitude; 
z/d 48, R 8500 
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Figure 10: Radial Distribution of Disturbance Amplitude; 
z/a 56, R 5700, R 8500 


47 


ga7 


OUY 


vu3 


0uz 


We 


Figure 





+ 
0 
a 
+ 
ays + 
+ 
Hy 
+ 
& 
Q 
+ BO HZ 1/2.1N Amp 
= 1) 20 H2 4/2 1N AMP 
ENG 
Nap ttt 
OSC. > as 
| 
926 033 oc8 00g 0!) 
V/a 
11: Radial Distribution of Disturbance Amplitude; 


z/ad 16, R 5700 


48 


O67 
+ 


a 
- 
é 4 J & + 
be ter 
8 + 
§ - 
U . i 
B , + 
0 als qj 4 * 
2 0 
+ - 
ue 
3 | + 30 HZ 4/2 IN AmP 
D 20 HZ 4/2 IN AMP 
Cj 0 
a 0 
° ooo o oponoooo -« 
5 | 
q 
S| 
008 053 OC6 cog Gia 


Y/o 


Figure 12: Radial Distribution Of Disturbance Amplitude: 
z/d 48, R 5700 
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Figure 13: Radial Distribution of Disturbance Amplitude; 
z/ad 16, R 8500 
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Figure 14: Radial Distribution of Disturbance Amplitude; 
z/a 48, R 8500 
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Figure 15: Radial Distribution of Disturbance amplitude: 
z/d 16, R 5700 
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Figure 16; Radial Distribution of Disturbance Amplit d 
z/a 32, R 5700 sheer 
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Figure 17: Radial Distribution of Disturbance Amplitude: 
z/d 16, R 5700 
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Figure 18: Radial Distribution of Disturbance Amplitude; 
z/a 48, R 5700 
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Figure 19: Radial Distribution of Disturbance Amplitude; 
z/ad 32, R 8500 
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Figure 20: Radial Distribution of Disturbance Amplitude; 
z/ad 48, R 8500 
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Figure 21: Radial Distribution of Disturbance Amplitude; 
z/a 32,R 8500 
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Figure 22: Polynomial Approximation to the Velocity Profile; 
z/ad 16, R 5700 
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Figure 23; Polynomial Approximation to the Velocity Profile: 
z/a 56, R 5700 
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Figure 2h: Polynomial Approximation to the Velocity Profile; 
z/a 16, R 8500 
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Figure 25: Polynomial Approximation to the Velocity Profile; 


z/a 56, R 8500 
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Figure 26: Decay of Disturbance Amplitude, R 7600 
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Figure 27: Neutral Stability Curve 
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